On Zermelo'-like problems: a Gauss- Bonnet inequality and 

a E. Hopf theorem 
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I 

^ ' The goal of this paper is to describe Zermelo's navigation problem on Ricman- 

nian manifolds as a time-optimal control problem and give an efficient method in 
order to evaluate its control curvature. We will show that up to change the Rie- 
mannian metric on the manifold the control curvature of Zermelo's problem has a 
simple to handle expression which naturally leads to a generalization of the classical 
, Gauss-Bonnet formula in an inequality. This Gauss-Bonnet inequality enables to 

' generalize for Zermelo's problems the E. Hopf theorem on flatness of Riemannian 

, tori without conjugate points. 
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1 Introduction 

> 

\^ I In the present paper we study a special class of time-optimal control problems on two- 

■ dimensional manifolds: the Zermelo'-like problems. By Zermelo'-like problems we mean 

the class of time-optimal control problems formed by the classical Zermelo's navigation 
I problems on Riemannian manifolds and the corresponding co-problems. 

' Our first goal in this paper is to describe these two problems and give an explicit 

. expression for their control curvature, the latter being the control analogue of the 

Gaussian curvature of surfaces. This is the purpose of Section [3l 

Zermelo's navigation problem aims to find the minimum time trajectories in a Rie- 
mannian manifold (M, g) under the influence of a drift represented by a vector field X. 
I The study of Zermelo's navigation problem began in 1931 with the work by E. Zermelo 

|22j and a while latter by C. Caratheodory in [T^]. In a recent paper [TU], Zermelo's 
navigation problem has been studied has a special case of Finslerian metrics and has 
been an efficient tool in order to give a complete classification of strongly convex Ran- 
ders metrics of constant flag curvature, the latter being the Finslerian analogue of the 
Riemannian sectional curvature. 

The co-Zermelo's navigation problem on a Riemannian surface {M,g) with drift T, 
where T is a one-form on M, is a time-optimal for which the maximized Hamiltonian 
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function h resulting of Pontryagin Maximum Principle has level one equal to the fiber 
bundle Ug(=A/{T"(j + Sq*}, where Sq* is the unitary Riemannian cosphere of the metric 
g. What is surprising with this definition of the problem is that it naturally leads 
to choose a good system of coordinates in which the control curvature has very nice 
and simple expression as a function of the drift one-form and the Gaussian curvature 
of the metric g. Contrary to the Zermelo's navigation problem the curvature of the 
co-Zermelo problem is much more readable than the one Zermelo problem itself and 
thus, much more exploitable. 

Another surprising property of the co-Zermelo's problem is that its flow is just a 
time rescaling of the magnetic flow of the pair {g,dT), the latter being solution of 
a fixed time variational problem. In particular, it implies that the curvature of the 
problem of a charged particle in a magnetic field is just a reparametrization of the 
curvature of co-Zermelo problem. 

We prove constructively that Zermelo's navigation problem on (M, g) with drift 
vector field X is feedback equivalent to a co-Zermelo problem but this time with respect 
to another Riemannian metric. This is the contents of Proposition [33] and its Corollary 
13.51 in ^3.4[ This proposition is fundamental because it points out that there are two 
different Riemannian metrics canonically associated to a given Zermelo problem. In 
particular, it implies that the two problems have the same curvature and also allows to 
see a Zermelo's navigation problem as its dual co-Zermelo problem and vice versa. 

This is of particular interest because the presentation of a given Zermelo's navigation 
problem as its feedback equivalent co-Zermelo problem has the serious advantage to 
present the curvature of the considered problem in an easier to handle formula. It is 
thus, clearly, the way to do. It also show how the classical Zermelo's navigation problem 
is linked to magnetic fiows. 

It is the second goal of this paper to show that there is a natural way to generalize the 
classical Gauss-Bonnet formula for Riemannian surfaces to an inequality for Zermelo- 
like problems. More precisely, we will see that, given a Zermelo'-like problem on a 
surface M there exists a canonically defined positive function (f> such that J"^ (pK dC ^ 
x(M), where x{^) is the Euler characteristic of M, Ti is the hypersurface and 
dC is the Liouville volume on 7i. Moreover, the function cj) takes the constant value 
equal to one if and only if the Zermelo problem is indeed Riemannian, in which case 
the inequality turns to be the classical Gauss-Bonnet formula. This Theorem 14.21 in ^ 
of the paper. 

Our last goal in this paper is to generalize for Zermelo's problems the E. Hopf's 
theorem which asserts that two-dimensional Riemannian tori without conjugate points 
are flat. This will be done in two steps following the Hopf's method. First we show 
that if a control system on a compact surface without boundary has no conjugate 
points then, its total curvature k dC must be negative or zero and, in the latter 
case its curvature must be zero identically. This is Theorem 15. II The second step is to 
use the Gauss-Bonnet inequality together with Theorem 15.11 to deduce flatness. In the 
Riemannian situation Theorem 15. II together with the Gauss-Bonnet inequality (which, 
in this case, reduces to the classical Gauss-Bonnet formula) imply straightforwardly 
flatness for tori without conjugate points. Of course Theorem 15.11 applies to Zermelo'- 
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like problems but, due to the presence of the function (p in the Gauss-Bonnet inequality, 
the situation is more delicate and essentially different. Indeed, a Zermelo'-like problem 
without conjugate points on a Riemannian torus is not necessarily flat unless its total 
curvature is zero. This situation is described in Theorem 15.21 and its Corollaries 15.31 
and 15.41 

To conclude our paper we discuss the further generalizations of the presented results 
to more general situations than the Riemannian one. We will see that even in the special 
case of Landsberg surfaces surfaces not all results can be transposed. 

2 Curvature of two-dimensional smooth control systems 

In the present paper smooth objects are supposed to be of class Let us fix some 
notations. For a two-dimensional manifold M, n : T*M — > M is the cotangent bundle 
to M. We denote by s the canonical Liouville one-form on T*M, sa = Aott*, A G T*M. 

If M is endowed with a Riemannian structure g, {■,-)g and \ ■ \g denote the Rie- 
mannian scalar product and the Riemannian norm respectively. Since the Riemannian 
structure defines a canonical identification between the tangent and cotangent bundle 
of M, we use the notations of the scalar product and norm indifferently for vectors and 
covectors, vector fields and one- forms. We denote by and S^* the unitary spherical 
bundle {v £ TM : \v\g = 1} and the unitary cospherical bundle G T*M : = 1} 
respectively. 

2.1 Definition 

We briefiy recall some facts concerning the curvature of smooth control systems in 
dimension two. For more details on the subject we refer the reader to one of the 
following items [5l[l9l[20]. 

Consider the following time-optimal smooth control problem 

q = f{q,u), qeM, u e U, 

q{0)=qo, q{ti)=qi, (2.1) 
ti min, 

where M and U are connected smooth manifolds of respective dimension two and one. 
For the above time-optimal control problem we denote hy h = max^g;/ (A, /(q, u)), 
A G T*M, q G M, the (normal) Hamiltonian function of PMP (Pontryagin Maximum 

Principle), hy Ti = h~^{l) C T*M, and by h the Hamiltonian field associated with 
the restriction of /i to Ti. Recall that the maximized Hamiltonian /i is a function on 
the cotangent bundle T*M one-homogeneous on fibers and non-negative. Under the 
regularity assumptions of strong convexity 

/(.,.)A^/0, ^A^!^/0, ,GM, uGf/, (2.2) 
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the curve 7iq = 7i f] T*M admits, up to sign and translation, a natural parameter 
providing us with a vector field Vg on Tig and by consequence with a vertical vector 
field V on TC. Vector field v is characterized by the fact that it is, up to sign, the unique 
vector field on TC such that 

LIs\'h = -s\n + bL^s\n, (2.3) 

where 6 is a smooth function on the level 7i. The function b, which is by definition a 
feedback- invariant, is called the centro-affine curvature. 

The vector fields h and v which are, by definition, feedback-invariant satisfy the 
following nontrivial commutator relation: 



V, h 



KV, (2.4) 



where the coefficient k is defined to be the control curvature or simply the curvature 
of the optimal control problem ()2.ip - ()2.2p . The control curvature is by definition a 
feedback-invariant of the control system and a function on TC (and not on M as the 
Gaussian one). Moreover, k is the Gaussian curvature if the control system defines a 
Riemannian geodesic problem. 

Example 2.1. Consider the time-optimal control problem corresponding to the geodesic 
problem on a two-dimensional Riemannian surface (M,g): 

q = u, q e M, M G TgM, \u\g = 1, 
q{0) = qo, q{ti) = qi 
ti — > min. 

In this case, the Hamiltonian function of PMP is given by 

hg{X) = A G T*M, 

and the vectors fields hg and Vg on /i^^(l) by 

d d 

hg = cos6'ei -|-sin6'e2 + (ci cos6' -|- C2 sin^) — , Vg = — , 

do od 

where (61,62) is a local (7-orthonormal frame whose structural constants ci, C2 are 
defined on M by 

[ei, 62] = ciei + 0262, ci, C2 G C°°{M), 
and 9 is the parameter on the fiber h~^{l) n TgM = Sg* defined by 

(A, ei(g)) = cos6', (A, 62(9)) = sin 0. 
The Gaussian curvature Kg of the surface (M, g) is evaluated as follows: 

Kg = -cj - cl+ Lej^C2 - Le^Cl- (2.5) 

Of course, for the Riemannian problem the curvature depends only on the base point 
q £ M as one can see from formula (|2.5|) but in general this is not the case: the control 
curvature depends also on the coordinate in the fiber Tig and thus is a function on the 
whole three-dimensional manifold TC. 
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2.2 Reparametrization 

In practice, it is sometimes easier for computations to not to consider the curvature 
itself but some reparametrization of it. We will thus see how the curvature k changes 
under a reparametrization of time. So let t = T{t) be a reparametrization of time. 
Under this reparametrization the ODE ^ = h{\) changes as follows: 

dXoT _dXdT _ dT-^ 

Thus, reparametrizing time just means to consider the field h in the form 

h = -, 

where ip € C°°{7i) is a positive function whose primitive along the trajectories of h is 
the time reparametrization function. Let v be is the vertical field defined hy v = ^/(pv. 
Vector fields h and v satisfy the following non trivial commutator relation: 



h, 



V, h 



(2.6) 



Indeed, denote for simplicity ip = ^JTp then, we have: 





V, h 














V, h 


+ V 










V, h 


+ 







V, h 



+ '4P'Lj:aI} v,h +tp'^Lf^'ip 



h, V 



+ ip'^Llipv + 



= kv + ^h, 

where .^i, G C°°{7i). From the previous computation one infers that the curvature 
and its reparametrization k satisfy the following relation: 



(2.7) 



We call the function k defined by the relation ()2.6p the ip-reparametrization of the 
curvature k. 



3 Zermelo'-like problems 

3.1 Zermelo's navigation problem 

In his article [2^ of 1931 Ernst Zermelo formulates the following problem: 
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"In an unbounded plane where the wind distribution is given by a vector field as a func- 
tion of position and time, a ship moves with constant velocity relative to the surrounding 
air mass. How must the ship be steered in order to come from a starting point to a 
given goal in the shortest time?" 

For our purpose we assume that we are working on a Riemannian surface in the 
presence of a stationary wind distribution that we call drift. Zermelo's navigation 
problem thus consists of finding the quickest path (in time) of a point on a Riemannian 
surface (M, g) in the presence of a stationary drift modeled by an autonomous vector 
field X S Vec M. This time-optimal control problem is read 

q = X{q) +u, qe M, u G TqM, \u\g = 1, (3.1) 
qiO) = qo, qiti) = qi (3.2) 
ti — > min, (3-3) 

and we call it Zermelo problem of the pair {g, X). The Hamiltonian function of PMP 
is 

h{X) = max ((A, X) + (A, n)) = (A, X{q)) + \X\g, (3.4) 

\u\g^l 

and the Hamiltonian vector field on TC = h~^{l) has the form 

h = X + hg+ (Umax, [ei,e2])g (nmax,^)^, + L[un..^^-vg] (^imax,^)g)f(,, (3.5) 

where the function Umax = ^^max(A) is the restriction to 7i of the maximized control 
obtained in the maximization ()3.4p . Relation ()3.5p leads naturally to an expression 
of the curvature of Zermelo's navigation problem p.ip - ()3.3p as a function of the drift 
X and the Gaussian curvature of the surface {M,g). We do not give here a precise 
formula for this expression of the curvature since it leads to a formula which is rather 
complicated and hardly exploitable except for very simple cases. We refer the reader 
to [2^ for a detailed description and coordinate expression of the curvature of this 
problem. 

3.2 Co-Zermelo's navigation problem 

Roughly speaking, whereas Zermelo's navigation problem was defined by its dynam- 
ics, i.e., as a subbundle of the tangent bundle over the state space M (in this case 
U^eM{^(g) C TM), co-Zermelo's navigation problem will be defined as a sub- 

bundle of T*M. Precisely, 

Definition 3.1. We call co-Zermelo problem of the pair (5, T) the minimum time 
problem for which the Hamiltonian function of PMP has level one equals to UggA/{Tg + 
cSf*} CT*M where T is a one- form on M such that \T\g < 1. 

Let h be the maximized Hamiltonian function of PMP associated to the co-Zermelo 
problem of the pair {g,T) which, let us recall it, is one-homogeneous on fibers and 
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non-negative. Denote by Ti the hypersurface h ^(1). By definition of the co-Zermelo 
problem the hypersurface 7i is characterized by 

(A-T,(;,),A-T,(;,)>^ = 1, VAGW. (3.6) 

Suppose now that A G T*M is a non zero covector such that h[\) ^ 0. Then, using the 
homogeneity of h we get 

Consequently, the covector A//i(A) has to satisfy equation (j3.6p . Plugging this covector 
in equation (13. 6p leads to 

(A-MA)T,(A),A-/i(A)T,(A)>' = MA)', A G r*M, (3.7) 

which gives an implicit definition for the Hamiltonian function h. Solving equation 
(I3T1) for h{\) gives 



,,,, -(^'T^^(A)>, + V(^'Tf'.(A)>' + (l-|T.(A)|^)|A|2 

^(-^) = ^ i-r 12 ' (3.8) 

where we have excluded the non-positive solution. We now derive the equation of the 
Hamiltonian field associated to h on the level surface 7i. If (p, q) is a canonical system 
of local coordinates on T*M and is the coordinate on fibers 7i H T*M, the vector 
field h is given by 



- - a a/i 

= TT^h + c— = —- 
oO opi 



d dh 
^ dqi dp2 



_d_ d_ 
n dq2 99 ' 



We define a parameter 9 on fibers 7i H T*^y^M in the following manner. Let (ei, 62) be 
a local g-orthonormal frame on M. Notice that equation p.6p can be locally rewritten 



(A - T,(A),ei(7r(A)))' + (A - T,(;,), e2(vr(A)))' = 1, VA G H. 
Hence, the fiber 7iq can be naturally parametrized by an angle 9: 

(A - T^(A),ei(7r(A))) = cos6', (A - T^(a), e2(7r(A))) = sin6'. 

In order to get the equations of the Hamiltonian vector field h, we write equation (j3.7p 
in coordinates (p, q) 

{P - HP: Q)'^q: f^l{<l)f + (P " KP: I)'^ qi ^2(9))^ = h^{P, q), 

and we differentiate it with respect to the pj's. We get 

^ ' dp dh \ dh 



k=l 
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Consequently, on the surface Ti. 

(p-T,ei)ei + (p-T, 62)6^2 



dpi 



^ l + {p-T, ei) (T, ei) + {p-T, 62) (T, 62) 
cos 9 e\ + sin 9 63 



1,2. 



l + cos0(T,ei) +sine(T,e2)' 
Thus the horizontal part of the field /i on H is 

- 1 / d d \ 

= — - ( (cos 9 e\ + s\n9 e\)— h (cos 6* e? + sin 6* 69)7: — ) 



— (cos 9 ei + sin 9 62) 



where 



93 (6l,g) = l + cos6l(Tg,ei(g)) + sin6' (T^, 62(5)) . 



Because h is the Hamiltonian field in restriction to 7i, we have ds\')-i{h, ■) = from 
which we can deduce the ^ of h. Let (e*,e2) be the coframe dual to (61,62) and 
denote sI-h = uj. In coordinates A = {9, q) on TC the Liouville one-form to takes the form 

io = (A, 61) el + (A, 62) 62 

= ((A - T, 61) + (T, 61)) el + ((A - T, 62) + (T, 62)) e* 

= cosset + sin 062 + T, (3.9) 

so that its exterior derivative is 

duo = — sin 9d9 A + cos 9d9 A 63 + cos 9de\ + sin 9de2 + dT . 

Using Cartan's formula for one forms di{X, Y) = Lx {C, Y) - Ly (C, X) - {C, [X, Y]), 
one easily see that 

del — ~cidVg, de2 = —C2dVg, 

where, as in Section [2?H ci, C2, are the structural constants of the frame (61,62) and, 
dVg = ej A 62 denotes the Riemannian volume element on M. Let Q € C°°(M) be the 
function defined by dT = —^IdVg and denote Cg = c\ cos 9 + C2 sm9. Summing up, we 
have 

duj = - sin 9d9 A el + cos 9d9 A 63 - {cg + n)dVg, (3.10) 

from which we get 



= du(h, •) = — csin 9el + ccos 060 + ^- — sin06i — ^^^-^ — cos 060. 



Hence, 

Cn + fl 

c = — . 
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Summing up, the Hamiltonian of the co-Zermelo problem reads 

1 / „ _ a 



h{e,q) = —r(f^ \ coseei{q)+s\nee2{q) + {cg{q) + ^{q)) 
or, equivalently 

/I = ^ fe% + , ^J(A) = 1 + (A-T,(A),T,(A)> , Ae?^, (3.11) 

where hg and Vg are defined as in Section 12.11 and is the diffeomorphism 

A ^ A + T^(A). 
Notice that (F^)-i = P-^. 

Remark 3.2. To conclude this section let us give a (coordinate free) formulation for 
the co-Zermelo problem of the pair (g, T) as a time-optimal control problem. According 
to (|3.11|) . this time-optimal control problem reads 

^ -, qeM, ue TqM, \u\g = 1, (3.13) 



(3.12) 



g(0) = qo, qih) = qi 
ti — > min, 

and the reader can check that the result of the maximality condition of PMP, max|u| =i (A, 
is the Hamiltonian function given by relation (j3.8p . 



3.3 Curvature of the co-Zermelo problem 

In order to get the expression of the curvature of the co-Zermelo problem, we first of 
all need to find the expression of the vertical field that satisfies relation p.3p . 
According to (f3l9]) and (fSTTO]) . 

to A dio = to = if'^g el A e*2 = ip^g dVy ^ 0, (3.14) 

which shows that (w, ^) forms a frame of horizontal one- forms on TC. The decompo- 
sition of the second derivative in this frame reads 



-—io + 



from which we deduce that the vertical vector field v that satisfies ()2.3p has the coor- 
dinate expression 

We now compute the curvature of the co-Zermelo problem according to relation ()2.4p . 
We find that 
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Proposition 3.3. The curvature of the co-Zermelo problem of the pair {g, T) is 



rr ifi I o TT 



r\2 
9 > 



(3.16) 



Proof. According to (IXTTD and (IXTHIl . 



h 



V 



h 



h = F^{hg + nVg), 



which imphes that it is enough for this problem to compute the (/jj-reparametrized 
curvature (defined in Section [2.2p . We have 



F7 



V, h 



9' 



Vq,h„ + nv„ 



ha 



Vg,hg 



Vg,hg 



[hg,Vg\ 9 



KgVg-nhg+Fj^^^^^.^Vg 



= {Kg + + L^^^^^^^n)F--v - nF-^h. 

According to ()2.7p the resuh fohows. 



We refer the reader to [18] for a detailed presentation of the co-Zermelo problem 
with linear drift on the Euclidean plane M^. In particular, using the reparametrized 
curvature, the author studied in great details the occurrence of conjugate points. 



3.4 Duality between Zermelo and co-Zermelo problems 

In this section we prove a proposition which asserts the feedback equivalence between 
the Zermelo and the co-Zermelo navigation problems. Although this proposition is 
simple indeed, it will have a fundamental role in the sequel due to fact that the curvature 
is much simpler to handle for the co-Zermelo problem than for the Zermelo navigation 
problem itself. 

Let {M,g) be a Riemannian manifold and fix an ^(-orthonormal frame (61,62). 
li X £ VecM, we define the local orthonormal frame for g associated to the vector 
field X with respect to the frame (61, 62) by 

6f = cos 61 + sin 62 



sin0'^6i 



where O-'^ (q) is the angle defined by 
'e^{q)=0 



cos 6*^(9) 



(X(g), 61(g)), 
\X{q)\g 



sm9^{q) 



cos ^^^62 



(X(g), 62(g)), 



1^(^)1^ 



ifX(g)=0, 
ifX(g)/0. 



(3.17) 
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In the same way if T G A^(M) we define the 5-orthonormal frame associated to the 
one-form T with respect to the frame (61,62) by 



-2 

3T 



eJ = cos ^^61 + sin ^^62 
62 = — sin 0^61 + cos ^^62 , 



where q 1-^ ^^('z) is the angle defined by 

-{q)=0 ifT, = 0, 

cos6l (g) = — , sm6l (g) = — if / 0. 



I 'v^ I 1 \M / I 



Notice that in this frames 



X = (x,6f )^6f = |xi,6f , T = (T,6r)6r = |T|,6r. 

Suppose for now that the Riemannian norm of the drift in our Zermelo navigation is 
strictly smaller than one. 

Proposition 3.4. Let (M, g) be a Riemannian surface. Let X be a vector field on 
M (respectively, T a one- form on M). There exists on M a new Riemannian metric 
9 = 9{9i^) (respectively g = g{g,T)) and a one-form T (respectively, a vector field 
X) such that the Zermelo problem of the pair {g,X) (respectively, the co-Zermelo 
problem of the pair [g, T) ) and the co-Zermelo problem of the pair (g, T) (respectively, 
the Zermelo problem of the pair {g, X) ) have the same Hamiltonians. 

Proof. Consider Zermelo's navigation problem (I3.ip - (l3.3p and let (61,62) be an or- 
thonormal frame for the metric g. Define some polar coordinates (p, 0) on the fiber 
T;M by 

p=\X\g, (A, 61) = p COS 6*, (A, 62) = /O sin 61, 
so that the Hamiltonian (|3.4|) takes the form 

h{p, 0,q)=p {\X{q)\g cos{9 - ^^(g)) + 1) , 

where 9^{q) is the angle defined by ()3.17p . Thus, the curve Tiq = /i~^(l) n T*M has 
the polar equation 

'^'^ = |X(g)|,cos(/-^-(g)) + l- ^^'^^^ 

Since \X\g < 1, the curve Tig is an ellipse centered at a focus. Moreover, this ellipse 
has for g a focal distance c = (p(7r + 6*^) — p{6^))/2 = \X\g{l — jXlp"-*^, a semimajor 

distance a = p{6-'^)+p{'K-\-6'^) = (1— |X|^)~-^, and a semiminor distance 6 = \J (P' — (? = 
(1-|X|2)-V2. 

In order to transform Zermelo navigation problem in a co-Zermelo problem, we con- 
sider the curve Hq as the drifted Riemannian cosphere at point q for a new Riemannian 
structure g on the manifold. In other words, we ask the one-forms 
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to form an orthonormal coframe for the new Riemannian structure g on the manifold 
and the one-form 



X 



to be the drift one- form of the co-Zermelo problem on (M, g). The corresponding (new) 
orthonormal frame (61,62) is characterized by 



((el, e^), (61,62)) = Id, 



which leads to 



ei 



[l-\X\l)ef, e, = Jl-\X\le-. 



Notice that we have (61,62) = (6-^^,62^) which shows in particular that \X\g = \T\g. 
The situation discribed above is illustrated by the picture below. 



'62(9) 





X{q)+SSrT^M 



Hg^Tg + Si* C T*M 



In order to complete the proof it remains to check that the Hamiltonian function 



(9.X) 



of the Zermelo problem of the pair {g,X) and the Hamiltonian function h. 



(9.T) 



of the co-Zermelo problem of the pair (g, T) are the same. For simplicity we denote 
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c = \T\g = \X\g. We have 

h^i'^\X) = (A, X) + |A|, = (A, cef ) + ^ (A, eff + (A, e^f 




(A, cei) + ^(A,ei)2 + (l-c2)(A,e2)' 
1^? 

(A, cei) + ^ (A, ei)2 + (A, e^f - (A, es)' - £2 (A, eif + c2 (A, ei)^ 

r^a2 

- (A, -cei) + ^( (A, ei)^ + (A, 62)^ ) (1 - c2) + ( - c (A, ei) f 

1 -C2 

-(A, (t, ei)et)g + ^|A|,- (1 - c^) + ((t, ei) (A, ei) f 

1-C2 

-(A,t),- + ^(l-|t|2-)|A|,- + (A,t)2- 
1-|T|| 

In order to prove the converse, one has just to permute the roles of vector fields and 
one forms in the previous considerations. ■ 

Zermelo's navigation problem and co-Zermelo's navigation problem which have the 
same Hamiltonian are said to be dual problems. The above proposition implies in 
particular that the two dual problems have the same curvature. This proposition can 
be reformulated as follows. 

Corollary 3.5. Two dual Zermelo's problems are feedback equivalent. 

Proof. Notations are these of the proof of the previous proposition. A similar com- 
putation computation as the one made in the previous proof shows that the two dual 
Zermelo's problems have the same sets of admissible velocities, i.e., that for every 
q £ M, {X{q) + u : u £ S^} = {u{l - {fg,u))-^ : u £ S^} (refer to equations (f3J]l 
and (j3.13p for the dynamics of Zermelo's problems). Thus, the feedback transformation 
u n(l — {lLq,u))~^ — X{q) has obviously the required properties. ■ 

3.5 Classical particle in a magnetic field on a Riemannian surface 

The motion of a charged particle of unit mass under the presence of a magnetic field 
is modeled by what is called the magnetic fiow. We will see here how the problem of a 
charged particle in a magnetic field is linked to the dual to Zermelo problem. Magnetic 
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flows were first considered by Arnold in [8] and by Anosov and Sinai in [7] but, it is 
Sternberg in [21] gave the first formulation of this problem using symplectic geometry. 

Let {M,g) be a two-dimensional Riemannian manifold and B S A^(M) a closed 
two-form thought as a magnetic field in which we have absorbed the electric charge of 
the particle as a parameter. 

The magnetic flow of the pair (^f, B) is the flow of the Hamiltonian 

with respect to the symplectic form as = ds + tt*B (see [21]). In the case where B 
derives from a magnetic potentiel, i.e., when B = dT, T S A^(M), the magnetic flow is 
also Hamiltonian with respect to the canonical symplectic form ds but this time with 
the Hamiltonian function 

1 1 

/imag(A) = - (A - T^(^), A - '^Tr{X))g = '^^gi'^ ~ ^7r(A))- 

A straightforward computation shows that the Hamiltonian vector field /imag associated 
to /imag in restriction to /imig(l) is given by 

^mag = F^{hg + ^Vg), 

where E C°°{M) is defined in same way as the function of the co-Zermelo problem. 
This shows that the equations of motion of a particle in a magnetic field are in fact 
the equations of motion of the reparametrized co-Zermelo problem. For this reason we 
define the curvature Kmag^' of the magnetic flow to be the c^J-reparametrized curvature 
of the co-Zermelo problem, i.e., 

K\^tr = ^g + ^' + L^n„.,f^ (3.19) 

so that, 

n'fJ' = {^ir' {i^' - 5(^J)) . (3.20) 

Remark 3.6. There is a theory on the reduction of the curvature of Hamiltonian flows 
by first integrals, see [3]. The reader can check that, what we have defined to be the 
curvature of the magnetic fiow corresponds to the reduced curvature of the Hamiltonian 
^mag on the level /imig(l)- 

4 A Gauss-Bonnet inequality for Zermelo's problems 

This section is dedicated to some global "Gauss-Bonnet properties" of Zermelo's prob- 
lems; key ingredients to prove Hopf's theorem for Zermelo problems (purpose of the 
next section). 

On the three-dimensional surface Ti there exists a canonical volume element, called 
Liouville volume element, defined by dC = — f\ds\'}-c. Since the Liouville one- form 
sIt^ is invariant by h so is dC, i.e.. 



Lr dC = 0, 



(4.1) 
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In the case of a Riemannian surface {M,g) the Liouville volume element on /i~^(l) = 
S*^ is called Riemannian volume element and we denote it by dTZg. In this particular 
case it is easy to check that dTZg is invariant by the vertical field Vg (actually the 
Riemannian case toghether with the Lorentzian are the unique ones satisfying the 
regularity assumptions ()2.2p for which the canonical vector field v, defined by relation 
(12. 3p leaves invariant the Liouville volume). Thus, being invariant by hg and Vg the 
Riemannian volume element is also invariant by their bracket, that is 

L[H„.,]dT^, = 0- (4-2) 

Using relation ()3.14p . one can easily checked that for the co-Zermelo problem of the 
pair {g, T) the two volume elements dC and dTZg are linked by the relation 

F^*dC = ^JoF-^dng, (4.3) 

where is the diffeomorphism defined by relation (|3.12|) . 

Lemma 4.1. Let {M,g) be a compact, orientable, two-dimensional Riemannian man- 
ifold without boundary. Let T be a smooth one-form on M. Then, 

^ j^^^ ° F'' d-Jlg ^ x{M), (4.4) 

j^^]t^':^TdC^x{M), (4.5) 

where xi^) Euler characteristic of the surface M . 

Proof. According to (|3.19|) . 

which, according to the classical Gauss-Bonnet formula and relation ()4.2p . is equivlent 
to 

A / ^mag"' ° F^ dTZg = 47t\(M) + 27r / dVg ^ 47r2x(M). (4.6) 
This proves relation (j4.4p According to relations (|3.2Up and (|4.ip . we have 



dC 



"^^^^'V^^O. (4.7) 



Relation (j4.5p follows from (j4.6p and (|4.7p . which completes the proof. 
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Theorem 4.2. Let M be a compact, orientahle, two-dimensional Riemannian manifold 
without boundary. If n is the curvature of a Zermelo'-like problem then, there exists a 
canonically defined positive function (p which is identically equal to one if and only if 
the problem is Riemannian such that 



Proof. It follows straightforwardly from the previous lemma and Proposition 13.41 ■ 

It immediately follows from the above theorem that 

Theorem 4.3. Zermelo's problems having non positive not identically zero curvature 
do not exist on two-dimensional tori. 

Proof. We prove the result by contradiction. Let k be the curvature of a Zermelo'-like 
problem on a two-dimensional Riemannian torus and let (p be the function of Theorem 
(14. 2p . Suppose that k ^ 0. Since n does not vanish identically, there exists a point 
X GTi. such that k(A) < 0, which, in addition with the fact that (p is a strictly positive 
function implies that 4>n dC < 0. But this contradicts the Gauss-Bonnet inequality 
of Theorem ()4.2p which, in this case reads (pndC ^ in'^xij^) = 0. ■ 

Remark 4.4. Although the previous theorem is an immediate consequence of inequal- 
ity ()4.8p . we want to point out that this theorem also follows from a more general 
fact if "non-positive" is replaced by "negative" in its formulation. Indeed, the flow 
generated by the Hamiltonian of a smooth control system having negative curvature 
is Anosov (see Moreover, in the appendix to the paper by Anosov and Sinai [7|, 
Margulis proved that if an Anosov flow operates on a three-dimensional manifold then, 
its fundamental group has exponential growth. Therefore, an Anosov flow cannot be 
carried by a three-dimensional torus since the fundamental group of the latter is the 
free abelian group Z,^ which is known to have polynomial and not exponential growth 
(see e.g. [E]). Finally, one easily check that the hypersurface W of a Zermelo'-like 
problem (of course, whose drift has Riemannian norm strictly smaller that one) over a 
two-dimensional torus is diffeomorphic to a three-dimensional torus. 

It's not worth mentioning that the Gauss-Bonnet (j4.8p inequality becomes an equal- 
ity NOT ONLY if the problem is Riemannian. Indeed, 

Proposition 4.5. The Gauss-Bonnet inequality of Theorem \4-^ is an equality if and 
only if the drift is identically zero or the Gaussian curvature of the manifold is zero 
and the drift has constant Riemannian norm. 





Moreover, when cj) is identically equal to one relation is the classical Gauss-Bonnet 
formula. 
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Proof. It follows from Proposition 13.41 that it is enough to prove the result for the 
co-Zermelo problem of the pair (5, T). Let M = UqOq where the O^s are domains 
of local ^(-orthonormal frames and let (61,62) be such a frame. From relation ()4.6p we 
know that 

/ iplKil'rdC = A7T\{M) + 27r [ n'dVg+ [ ( ^-{' ) F~-'*dng (4.9) 
Jn Jm Jn \ ^V'g / 

so that the Gauss-Bonnet inequality becomes an equality if and only if 

^ = ^^^d L^^^^y,=0 (4.10) 

identically. On the one hand, the condition = means that the drift form T is closed 
(recall that 0, was defined by dT = QdVg), which implies 

= dT(ei, 62) = LeiT2 - Le2Ti - Tici - T2C2, (4.11) 

where Ti = (T,6i) and T2 = (T,62). 

On the other hand, keeping in mind that $7 = holds true, condition Lr ipj. = 

"■mag y 

reads Lp-^j^ (^J = 0. According to the notations of Example 12.11 that is 

= ^coseei+sm0e2+(cicose+C2sine)^(l + ^lCOSe + T2Sine) 

= {Le^Ti + C1T2) cos^ 9 + (-^62X2 - C2T1) sin^ 6 

+ (LeiT2 + Le2Ti - ciTi + C2T2) COS ^ sine (4.12) 

Equations (j4.1ip and (I4.12p are thus equivalent to the system of equations 

LeiT2 - LeaTi - CiTi - C2T2 = 

LejTi + C1T2 = 

Le2^2 - C2T1 = 

LejT2 + LeaTi - CiTi + C2T2 = 0. 

Replacing the first and last equations respectively by there sum and difference we 
equivalently get 



LeiT2-ciTi = (4.13) 

LeiTi + ciT2 = (4.14) 

ie2T2-C2Ti = (4.15) 

Le2Ti + C2T2 = 0. (4.16) 



Now we differentiate equation (|4.16p along 61 and subtract it the differentiation along 
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62 of equation ()4.14p . According to (12. 5p , we get 
= L^. Km - LeadUl]) 

= Lej o LgjTi + C2LeiT2 + T2LeiC2 — ° L^^Ti — CiL(,2^2 — T"2-^^e2Cl 

= -^^[61,62]^! + '^2{LexC2 - Le^Cl) + (TiCi)c2 - (TiC2)ci 
= CiLeiTi + C2Le2Ti + T2(LeiC2 - Le2Cl) 

= T2(-cf-ci + LeiC2-Le2Ci) (4.17) 

= T2K3. (4.18) 

In the same way, using this time equations (|4.13|) and (|4.15|) we get 

= Km - Lei (1^^ = TiKg. (4.19) 

If the Gaussian curvature is identicaUy equal to zero then, the Riemannian manifold 
is a flat torus. In this case we can chose local coordinates ((71,^2) on M such that 
ei = ^ and 62 = In these coordinates equations (|4.13p . (14.141) . ()4.15p and (I4.16P 
read 

Le,T, =0, i,j = l,2, 

which obviously implies that the coefficients Ti and T2 are constant. Therefore T has 
constant Riemannian norm. 

If the Gaussian curvature is not identically equal to zero then, it follows from equa- 
tions ()4.17p and (I4.19P that the form T must be zero wherever Kg is different from 
zero. Consider the set A = {q £ M : Kg{q) =0}. If the interior of A is empty it 
follows from its continuity that T vanishes identically on M. If the interior of A is non 
empty, a similar reasoning as above (done on successively on each domain Oa) and 
the continuity of T imply that T has constant Riemannian norm in restriction to the 
closure of the interior of A. But, T| jy/^^ = and since (Af \ j4) n clo int A 7^ 0, by conti- 
nuity we must have T = identically on M. This completes the proof of the theorem. ■ 



Remark 4.6. Equations (|4.10p are indeed equivalent to the unique equation Lr ipj; = 

"•mag y 

0. Namely, Lr (pi is a polynomial of degree two in cos 9, sin9. In particular we have 

"■mag y 

= L^ (^J(7r/2,g)+L^ ipl{-7r/2,q)=QT^ 

"-mag y ^ ' ''■mag ^ 

= L^ 99^(0, g) + L^ ^J(7r,g) = 17T2, 

"-mag y "-mag y 

which obviously implies that = 0. 



5 A E. Hopf theorem for control systems 

It is well known that Riemannian tori without conjugate points are flat. This theorem 
was first proved by E. Hopf in 1943 for the two-dimensional case (see [T7]) and for 
higher dimensional manifolds it was proved by D. Burago and S. Ivanov in 1994 (see 
|llj). We give in this section a generalization of Hopf's for control systems. 
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5.1 Jacobi curves 

We introduce here the Jacobi curves which are a generaUzation of the space of Jacobi 
fields along Riemannian geodesies. Since the construction of Jacobi curves does not 
depend on the dimension of the manifold, we begin with the general case to then go to 
our special low-dimensional case. 

Let h be the Hamiltonian function of PMP for a time-optimal smooth control prob- 
lem and H its hypersurface /i^^(l). Let e*'* : H —>^ H denote the flow generated by 
the Hamiltonian field of PMP h. This flow defines a one-dimensional foliation JF of 

— ^ 

7i whose leaves, the trajectories of h, are transverse to the fibers T*M, q £ M. This 
foliation enable us to make the following symplectic reduction. 
Consider the canonical projection 

7f : ?^ ^ s = n/r. 

The quotient space E, space of trajectories of h, is, at least locally, a well-defined 
smooth manifold and carries a structure of symplectic manifold with symplectic form 
a characterized by the property that its pull-back to TC is the restriction crly^. 

Let n C T7i denote the vertical distribution, i.e., IIa = ^^'^^^(a), X £7i. The curve 

Ja : M ^ r^(A)S 

is called Jacobi curve at A. Because the Hamiltonian flow preserves the symplectic 
structure, it is easy to check that the spaces Jx{t), t G M, are Lagrangian subspaces 
of the symplectic space r^(A)^ so that the Jacobi curves are curves in the Lagrangian 
Grassmannian L(r^(A)^)- 

Recall that the Lagrangian Grassmannian L(T^(a)^) of the symplectic space T^(a)S 
is defined by: 

L(r^(A)S) = {A C r^(A)S I = A}, A^ = G r^(A)S | A) = 0}. 

The Lagrangian Grassmannian of a symplectic space is a well-defined smooth and com- 
pact manifold. In our particular case of a two-dimensional manifold M, the Lagrangian 
Grassmannian L(r^(A)^) is diffeomorphic to the one-dimensional real projective space 
MP(1). Moreover, since the vertical distribution H is generated by the vertical vector 
field V the Jacobi curve can written as 

JaW = M(7f,e*"<^'^^;(A)) . (5.1) 
We say that a point e*'*(A) is conjugate to A (or time t is conjugate to zero) if 

Jx{t) n Ja(o) / {0}. 

Most of the material presented in this section can be fund in great details in the 
papers [H [H [6]. 
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5.2 A E. Hopf theorem 

In this section we prove the following 

Theorem 5.1. Consider a control system q = f{q,u) on a compact surface M without 
boundary. Assume that the curves of admissible velocities are strongly convex curves 
surrounding the origin. Then, if there is no conjugate points on M the total curvature 
Jy^ K dC must be negative or zero. In the latter case k, must be zero. 

Proof. Notice that because the curves of admissible velocities are strongly convex 
curves surrounding the origin, the manifold TC is compact. Although the proof we 
make here essentially follows the one given by Hopf in [T7], it will however be exposed 
in a more intrinsic and geometrical manner. The first step in the proof consists in the 
construction of a well-defined function on any extremal of our system, i.e., a function 
that does not depend on time but only on the point of the extremal. To do so we use 
the notion of Jacobi curve described in the previous section. 

Let A be a point of the hypersurface 7i C T*M and let J\{t) be the Jacobi curve 
associated with the extremal e*'*(A). we have 

Jx{t) = R (Tt^e'^^^^viX)^ £ MP(1), 

with 

e*^^'^^(A) = Pit, X)v{X) + 7(t, A) \v, h] (A). 

Considering (/3 : 7) as homogeneous coordinate in ]RP(1), we can identify the Jacobi 
curve with the curve 

t^m,X):j{t,X)). 

From the non existence of conjugate points it follows that 7(t, A) 7^ for t ^ 0. We 
can thus use the chart (/3 : 7) 1— > ^ and make the identification 

It turns out (see e.g. [SJ [19]) that the coefficients /3 and 7 are solutions of the Cauchy 
problems 

(3 + Kt(3 = 0, m = h /?(0)=0, Kt = K{e'^{X)), 
7 + Kt7 = 0, 7(0) = 0, 7(0) = 1, 

which shows in particular that (3 and 7 are two linearly independent solutions of the 
Hill equation x + Ktx = 0. The derivative with respect to time of the function yt is 

dyt ^ $1 - Pi 
dt ~ 72 

and because the Wronskian 



/3(0,A)7(0,A)-/3(0,A)7(0,A) = -1, 
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the function yt is strictly decreasing or, equivalently the Jacobi curve is strictly de- 
creasing in MP(1). Since yt is strictly decreasing its limit as t goes to infinity exists. 
Moreover, because of the non existence of conjugate points, this limit is finite. Indeed, 
notice that because of the initial conditions /3(0, A) = 1, 7(0, A) = and 7(0, A) = 1 we 
have for t small enough 

yt{\) > 0, y_t(A) < 0. (5.2) 



So if we suppose that 



lim yt{\) 

i— >+oo 



-00, 



(5.3) 



it would follow from Equations (15. 2p and from the strict monotonicity of yt the existence 
oft" < < t+ such that yt-{^) = yf+(A). Then, the time reparametrization t = t — 
would imply that time t = — is conjugate to r = 0, which is a contradiction. 
Hence, the function ?/+ defined by 

y+(A) = lim yt{\), XgH 

I— > + oo 

is a well defined function on the manifold 7i. Equivalently, the distribution 11^ S T7i 
defined by 



lim Jx{t) 

t—t + OO 



y^v + 



V, h 



is a well defined distribution on 7i transverse to the vertical distribution. This distri- 
bution is, by definition, invariant by the fiow of h. In terms of function this 
invariance reads 



D, h 



a y V + 



V, h 



+ 



V, h 



h, y V + 

or, equivalently 

LfiV^v + y 

where a is function on 7i. Solving (j5.4p for a gives 

a = — y+ and Lj^y^ + k — ay 
which shows that y+ satisfies the Riccati equation 

+ + « = 0. 



V, h 



ay V + a 



0, 



V, h 



(5.4) 



(5.5) 



As a limit of smooth functions, y* is clearly measurable, y* is also uniformly bounded 
as shows lemma 2.1 of |14] and thus it is integrable. If we now integrate equation (j5.5p 
over 7i with respect to the Liouville volume dC, the first term in the left-hand side of 
(|5.5p will disappear since the Liouville volume is invariant by the flow of h. As a result 
we obtain 



K dC 



n 



y+^dC 



H 



(5.6) 



which immediately proves the validity of the first part of the theorem. If we now sup- 
pose that the total curvature k dC is zero it follows from (j5.6p that the function y+ 



U. SERRES 



22 



must vanish everywhere on 7i. According to (I5.5p n must therefore vanish everywhere. 



We say that a control system q = f{q,u) is flat if it is feedback equivalent to a 
control system of the form q = f{u). 

In the Riemannian case, a direct consequence of the Gauss-Bonnet and Theorem 15. II 
is that two-dimensional Riemannian tori without conjugate points are flat. Contrary to 
the Riemannian situation, we shall see that Zermelo'-like problems without conjugate 
points on tori are not necessarily flat. 

The following three theorems give us a well understanding of the Zermelo'-like 
situation. To simplify notations, we omit the pair {g, T) in the writing of curvature 
and, the diffeomorphism ()3.12p in formulas since, anyway, its action is clear. 

Theorem 5.2. Consider a co-Zermelo problem on a compact Riemannian surface with- 
out boundary. If there is no conjugate points then the total curvatures n^azdC and 
fgg* K-msLgdTZg have to be negative or zero. In the latter case the considered co-Zermelo 
problem is flat. 

Proof. The part of the theorem concerning k^^z is given by Theorem 15.11 In order to 
check that the same conclusion holds for the curvature Kmagj let us see how changes the 
function y+ constructed in the proof of Theorem 15.11 under a reparametrization. For 
simplicity, denote = In a general manner, we have 



and we compute the new function y+: 

y'^v-\- v,h = y+'i/'-O + ipv ^ 

L J L . 



and V = ipv, 



y+ijv + \ ^ 

i' 



V, h 



V, h 
mod h 



mod h 



We thus have 



y V + 



V, h 



y+ = y+ij)"^ - V'L^V- 



(5.7) 



In the same way as for the function ?/+ it is easy to see that the function y^ satisfies 
the Riccati equation 



Notice that the Riemannian volume element dTZg is invariant by Aimag since 



(5.^ 



dUg = L^vgdlZg = d{ndV) = d{dT) = 0. 



Therefore the integration of ()5.8p leads to 



59* 



l^msugdl^g 



y+^dlZg < 0. 



53* 
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This prove the first part of the theorem and a similar argument as the one used in the 
proof of Theorem 15.11 shows that 

'^mag is zero GVGrywtlGrG wtlGn J^^* ^mag 

dUg = 0. 

We now complete the proof showing that the co-Zermelo problem is flat when the 
total curvatures k^oz dC and /^g, ^mag dTZg are both zero. In that case, we must 
have Kmag = and k^^z = everywhere. In particular it implies 

= / Kmag dTZg = I (fK.^z dC = 0, 

Jss* Jn 
i.e. (see the proof of Lemma l4.ip . 

= 47r2x(M) + 2^ / n^dVg = 4Tr\{M) + 2Tr [ n^dVg+[ ( ^] —, 
JM JM Jn \ / V 

which is equivalent to Lj^if = 0. Therefore, according to Theorem 14.51 and Remark 14. 6^ 
or the form T is different from zero and in this case the conclusion is obtained, or the 
form T is identically zero and in this case the problem is Riemannian. In the latter case, 
we have = Kmag = K-g which, on the one hand, implies that the Riemannian surface 
is flat and, on the other hand, according to the Gauss-Bonnet formula, it implies that 
the surface is a torus. The proof is complete. ■ 
A direct consequence of this theorem are the following corollaries. 

Corollary 5.3. // a co-Zermelo problem on a two-dimensional Riemannian torus has 
no conjugate points then, the torus is flat and the drift one-form is closed. In particular, 
time-optimal trajectories are straight lines. 

Corollary 5.4. Zermelo'-like problems without conjugate points on two-dimensional 
Riemannian tori are flat if and only if their total curvature is zero. 

5.3 A natural question 

In the proof of Theorem 15.11 we constructed a function y+ well-defined on TC that sat- 
isfies Riccati equation (j5.5p . This construction is valid along every regular extremal 
without conjugate points. Recall moreover that a control system with negative curva- 
ture does not admit conjugate points. A very natural question is thus the following: 

Does a control system without conjugate points admits a non positive (p-reparametrized 
curvature ? 

Since the function y+ satisfies Riccati equation (|5.8|) . the question can be reformu- 
lated in the following manner: does there exists a non vanishing function ip, say ip > 
for simplicity, such that Lfjj^ = 0, or equivalently such that L^y'^ = 0? According to 
relation (15. 7p . 
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SO that (dividing by V'^) Lj^y'^ = is equivalent to 



i.e., to 



Lj;y+ + 2y+L^ log V - (L^ log V^)^ - = 0. (5.9) 



Denote (7 = logV'- We have 



or equivalently 

This implies that equation (|5.9|) is equivalent to 

Lgy+ + 2y+Lg5 - 2(L;jg)2 - llg = 0, 
i.e., to the Riccati equation 

L^z + 2z^ - 2y+z - Lj^y+ = 0, (5.10) 

where we have set z = Lj^g. 

The function z = y+ is solution to Riccati equation (j5.10p . Thus we will have the 
required reparametrization of h if we can solve the equation 

4log^ = y+ (5.11) 

globally on the three-dimensional manifold TC. The first thing we need for the resolution 
of equation (|5.1ip is the continuity of the function y+ on TC. In the case of hyperbolic 
systems (see [16] for the definition), the function y+ is easily seen to be continuous 
due to some "exponential estimates" along the stable distribution (see [IS])- Also, for 
such systems the function is in general never differentiable and even never Lipschitz 
continuous but only Holder continuous (see [IH] Theorem 19.1.6 of Chapter 19). In the 
case of systems without conjugate points the situation is quite different because we do 
not have the exponential estimates and by consequence the continuity of the function 
is not so obvious. What we can ensure is the following. 

Lemma 5.5. The function y+ defined above is upper semi- continuous. 

Proof. Let (An)neN C W be a converging sequence to X £ 7i. Since yt{Xn) is decreasing 
in t, it follows that 

2/t(An) ^ 2/^(An) = lim yt{Xn)- 
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Taking the liminf as n tends to +00 in the previous relation, we get since ?/t(A) is 
continuous in (t, A) 

yt{\) ^ liminf y+, 

An — >X 

and then, letting t going to +00 leads to 

y^(A) ^ liminf , 

An — *A 

which proves the upper semi-continuity of y+. ■ 

Suppose that the function is continuous. It implies that we can solve locally 
equation ()5.1ip . In order to solve this equation globally, the question is more delicate 
because the problem is closely related to the fact that the quotient manifold S (defined 
in Section 15. ip is globally defined. It is not our scope to discuss this problem here. 
However we can say the following. Let M be the universal covering of M. because of 
the non existence of conjugate points, M is diffeomorphic to M?. Let 

h=f{q,u), q£M, ueU, (5.12) 

be the lift on M of the control system q = f{q, u), and 7i be the corresponding Hamil- 
tonian hypersurface. Then, the continuity of y+ implies that when the control system 
q = f{q,u) has no conjugate points then, there exists a reparametrization of h or, 
equivalently globally defined function ^ satisfying equation (jS.lip . such that the lifted 
system (I5.12p has negative curvature. Unfortunately, is in general not a continuous 
function as shown by Ballmann, Brin and Burns in [9] where they gave an example 
of a two-dimensional compact surface without conjugate points where y+ fails to be 
continuous. 

A nice work would be to characterize two-dimensional smooth control systems with- 
out conjugate points where this function fails to be continuous. What is the geometrical 
property that prevents y+ from being continuous? 

6 Conclusion 

We conclude this paper with a brief discussion of the extension of our results to more 
general structures than Riemannian surfaces. Of course, Zermelo'-like problems can be 
defined on any manifold equipped with a geometric structure defined by an optimal 
control problem of type (|2.ip - (|2.2p . A natural class of geometric structures on which 
generalize our result is the class of manifolds equipped with a Finsler metric (see the 
book of Chern and Shen [13] for a nice and brief presentation of Riemann-Finsler 
geometry). Unfortunately, since the Gauss-Bonnet formula is not true for any Finsler 
surfaces results from Section [5] can not be extented to all of these structures. One has 
to limits itself to Zermelo'-like problems on Landsberg surfaces on which almost all 
results from Section U] remain true. Roughly speaking, a Landsberg surface is a Finsler 
surface on which the Gauss-Bonnet formula remains true, (up to change the classical 
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27r factor in the formula by the centro-affine length (. of the curve 7iq which, let us recall 
it, is defined by £ = lAUq where /i is a one-form on the hypersuface 7i such that 
(/X, = 1 ). Without entering into details one can see that the Gauss-Bonnet formula 
still holds on Landsberg surfaces due to the fact that the centro-affine length of the 
curve 7iq does not depend on the base point q. This property is characterized by the 
fact that the invariant h that appears in relation <\2.2>\) is a first integral of the vector 
field h (see |19] for details). If we now consider Zermelo'-like problems on Landsberg 
surfaces, on the one hand, the Gauss-Bonnet inequality (14. 8p still holds true. The proof 
is the same but this time one has to be more carefull because the Landsberg volume 
element c?>Ciand is not invariant under the vertical Landsberg field I'land- Indeed, one 
can easily check that (i£iand = ^dC\s,^A- Anyway, (i£iand is still invariant under 

the bracket [^land) '^land] since 

L\r l^Aand = {Lr b)dC\a,ad = 0. 

["■landil'landj ^ "•land ' 

On the other hand, Theorem 15.21 and its Corollaries 15.31 and 15.41 do not generalize to 
Zermelo'-like problems on Landsberg surfaces. The reason is the following: Landsberg 
surfaces of zero curvature are not necessary flat (see [TH], Theorem 4.3.3). 
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